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Abstract 

An atom recoils when it undergoes spontaneous decay. In this paper we present a microscopic 
calculation of the recoil of a source atom imbedded in a dielectric medium. We find that the source 
atom recoils with the canonical photon momentum n/i&o, where n is the index of refraction and 
hko is the photon momentum calculated at the source atom atomic frequency loq. We also show 
explicitly how the energy is conserved with the photon inside the medium. 



I. INTRODUCTION 



The momentum of a photon in a dispersive medium has been considered by many authors, 
due to its conceptual and practical importance. One of the issues is how the momentum is 
conserved when a photon with momentum Kk in the vacuum is scattered from an atom in 
the medium to a new momentum state Kk' in the vacuum. Should the momentum imparted 
on the atom be the difference of the momenta in the vacuum Kk' — Kk or the difference of 
the canonical photon momenta in the medium n (Kk' — Kk) [6]? One might argueQ that, 
assuming the medium is dilute, the atom is localized in the vacuum space between particles 
of the medium, therefore the photon, before and after it hits the atom, travels in the vacuum 
and momentum should be conserved in terms of the vacuum momenta. On the other hand 
jlj], one can also argue that it is the macroscopic field of the incident wave that induces and 
interacts with polarizations and polarization currents of atoms in the medium and therefore 
the imparted momentum should be the difference of the canonical photon momenta in the 
medium. Experimentally, this issue has been studied in two systems. One measures the 
recoil of a mirror immersed in a liquid when the light is reflected from it |3j, and more 
recently j^J , a measurement of the recoil frequency of the Bose condensed 87 Rb using a two- 
pulse Ramsey interferometer. Both experiments confirm that atoms recoil according to 
canonical photon momenta. 

Most theoretical studies related to this issue deal mainly with classical fields. Some of 
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the recent work by Loudon 4] and Nelson p clarified some issues related to momentum in a 
dielectric from a quantum and microscopic perspective. Milonni and Boyd |6] consider a case 
where a source atom imbedded in the medium recoils due to its spontaneous decay. They 
find that the source atom recoils according to nhu>o/c, where uq is the atomic frequency and 
c is the speed of light in the vacuum. Their calculation is based on a Heisenberg Picture 
approach, where the operator expectation value, (P 2 ) is calculated to be r^Ktul/c 2 . Here we 
present a similar calculation in the Schrodinger Picture. The calculation in the Schrodinger 
Picture is particularly revealing because it includes explicitly processes that are responsible 
for the modification of the momentum imparted by the photon. As we show in the following, 
the photon travelling in the medium experiences a series of scatterings from medium atoms. 
Different scattering amplitudes interfere to shift the central frequency of photons in the 
field. Since the source atom is coupled directly to the field, by momentum conservation, the 
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source atom recoils according to the modified central frequency of photons. Our calculation 
is based on a quantum field quantized in the free space, which allows a separate description 
of the field and the medium, i.e. any wavelength and frequency of the field are calculated 
unambiguously in the vacuum. 

II. MODEL AND THE RECOIL CALCULATION 

The calculation is based on a model that we have used previously [?J • The source atom, 
with finite mass M, centered at position (R) = 0, has two internal levels, whose frequency 
separation is denoted by uj$. The uniformly distributed dielectric atoms have J = ground 
states and J = 1 excited states. The frequency separation of the ground and excited states 
is denoted by to. We assume the mass of the medium atoms to be infinite, which allows us to 
ignore recoil of the medium atoms. At t = 0, the source atom is excited to the m = excited 
state sublevel with center of mass momentum (P) = 0, and (AP 2 ) <C (hk ) 2 , the dielectric 
atoms are all in their ground states, and there are no photons in the field. The process we 
consider one in which radiation emitted by the source atom is scattered by dielectric atoms. 
The medium is modeled to be infinite, i.e. photons are always inside the medium. The 
vaccum field amplitudes, the medium atoms' excited state amplitudes, and the source atom 
center of mass motion is caculated as t —>■ oo. It is assumed that the medium atoms are far 
detuned from the source atom, uj u , and also uo ^> 7, where 7 is the spontaneous decay 
rate of the source atom. 

The free part of Hamiltonian describing such a system is 

#0 = ^y<t z + Yl ^Y a z j) ( m ) + ^ a ix a k\ + 7^ (!) 

j m=-l 

where a z = (|2) (2\ — |1) (1|), \2) and |1) are the m = excited and J = ground state kets of 
the source atom, respectively, <7z\m) = (jm)^ (m\ — \g)^ (g\j is the population difference 
operator between excited state \ J — l,m) and ground state | J — 0, g) of dielectric atom j, 
and OkA is the annihilation operator for a photon having momentum k and polarization A. 
We have also included a term describing the external motion of the source atom, where P 
is the momentum operator for the source atom. A summation convention is used, in which 
any repeated symbol on the right hand side of an equation is summed over, unless it also 
appears on the left-hand side of the equations. The interaction part that couples the field 
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with the atoms is 
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where <r± are raising and lowering operators for the source atom and (t± (m) are raising and 
lowering operators between the excited state \ J = l,m) and the ground state \ J = 0,g) of 
dielectric atom j, and is the matrix element of the dipole operator for the source atom 
and \x m is the matrix element of the medium atom between levels \ J = l,m) and \J = 0,g). 
We have taken the reduce matrix element, //, of dipole moments of the source and medium 
atoms to be equal to simplify our calculation. The polarization vectors are defined as 

e^. = cos 6*k cos 0k x + cos 6*k sin k y — sin 6 k z (5) 
e^' = - sin k x + cos k y- (6) 

The source atom interacts only with the z component of the vacuum field. Since the medium 
atoms are far detuned from the source atom, we include anti-rotating terms in the field- 
medium atoms' interaction. We have not included such terms for the interaction Hamiltonian 
between the source atom and the field because we have chosen the initial state to be the 
source atom excited with no photon in the field. 

Instead of writing amplitude equations for different states and then identifying terms 
corresponding to contributions from different processes [8(, we adopt a resolvent approach 
0, which allows us to write amplitudes directly from diagrammatic representations of the 
scattering processes. We want to calculate the recoil energy of the source atom, which 
includes the contribution from three amplitudes: the source recoiled atom with one photon 
in the field, the recoiled source atom with one medium atom excited, and the recoiled source 
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atom with both the field and a medium atom excited. These amplitudes are represented by, 

b k = (l,q;0;k|E/(oc)|2,O;0;O) (7) 
b mj = (l )q ;m i; 0|[/(oo)|2,0; <?;()) (8) 



Jmjkk' 



(l,q;m,;k,k'|tf(oo)|2,0; <?;()), (9) 



where the state of the whole system is labelled with source atom internal states, source atom 
wave vector, medium atoms internal states, and quantum field wave vectors. For example, 
|1, q; g; k) is a state of the source atom in the ground state |1), with a momentum ftq, the 
medium atoms in the ground state \g), and a photon with a wave vector k present. We use 
\rrij) to label the state with a medium atom at position Rj being excited to the |J = l,m) 
state. The photon polarizations are not written explicitly in the formulas. It should be 
understood, in the following perturbative calculations, any intermedia states' photon po- 
larzations are summed over in amplitudes, while the final states' photon polarizations are 
summed over in probabilities. 

The evolution operator can be expressed in terms of the retarded propagator as 

1 f 00 

U(t) = / dE exp (-iEt/H) G (E) (10) 

and the resolvent operator, G (E + zO), is defined as 

1 



G(E) 



E + iO-H 

where H is the total Hamiltonian, E is the energy of the system, and the infinite small 
imaginary iO prescription is used to yield retarded propagation. In order to carry out a 
perturbative calculation in orders of the interaction V, we recast the resolvant operator in 
the following exact form 

G(E) = G {0) (E) + G {0) (E)VG(E) (11) 
Here G^ (E) is the zeroth order propagator 

G<°> (E) 1 



E + iO — H Q 

The processes that we want to take into account are shown in Fig(P). Let's focus on the 
calculation of b k first, 

i r°° 

— lim / dEexp(-iEr/h)(l,q;g m ;k\G(E)\2,0;g m ;0), (12) 
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(a) (b) 

FIG. 1: (a) shows that the excited source atom |2) spontaneously decays by radiating a normalized 
photon with wave vector k and goes to the ground state |1) with external momentum —Kk. The 
thick wavy line corresponds to the normalized photon propagator and striaght solid lines correspond 
to the atom propagator, (b) shows the photon of wavevector k being scattered by media atoms. 
We include both the rotating and anti-rotating contributions in the intermedium states. Medium 
atoms can get excited to J = 1 sublevel m by absorbing a photon k and then go to the ground 
states by radiating photons. They can also get excited by emitting photons and then deexcited 
by absorbing other photons. In this diagram, the time propagation is from the left to the right. 
Averaging over the positions of the medium atoms results in conservation of the photon momenta 
before and after the scattering from medium atoms. 



which requires a calculation of the matrix element of the operator. According to equation 
(111)) , the resolvant operator can be expanded to first order in V as, 

(l,q;g m ;k\G^\2,0;g m ;0) 

= (1, q; g m , k| G<® (E) |1, q; g m ; k) (1, q; g m , k| V |2, 0; g m , 0) (2, 0; g m , 0| G (0) (E) |2, 0; g m , 0) 

= %^q,-k 0™ ■ e k) p Z 1 , TZ IT— ^W mi0 . (13) 

E — nu> + inn^i E + iO — nuf — rkj k 

We neglect the possibility that the photon can be scattered by medium atoms in this or- 
der. In the calculation, the source atom recoil momentum equals the inverse of the photon 
momentum, q = — k (conservation of momentum), resulting from evaluating the matrix el- 
ement (p = hq\ e lk R |p = 0). We have rearranged the order in the last line of the above 
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formula ([13]) so that the first part, 

describes the decay of the source atom in the medium. The imaginary part in the dominator 
ihn'-f is added to represent spontaneous decay [9 ! ]. The modification of 7 in the vacume to 
rry in the medium is related to the process that photons scattered by the media atoms are 
reabsorbed by the source atom, as shown in papers QjjJ. The second part of Equation ((IB)) , 

' W m , , (15) 



E + iO — hw r q — huik 



describes the propagation of the photon, where uj r q = hq 2 /2M is the recoil frequency associ- 
ated with the emission of a phone with wavevector q and W m; o is the transverse polarization 
tensor of the photon defined as 



'e^OeM =(l-jfc®*) , (10) 
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where two relationships 



and 



t4 ■ Cfe = (£fc • Hm) W m,0 (17) 

have been used. 

In order to carry the calculation of the matrix element of the resolvant operator, 
(1, q; g m ; k| G \2, 0; g m ; 0) , to higher order, it is neccessary to consider processes that the 
photon is scattered by medium atoms. Including these processes modifies the photon prop- 
agator (^3). For a dilute medium with iVAg <C 1, where N is the density of the medium and 
Ao = 2ttc/ujo is the photon wavelength, it is appropriate to make use of the independent 
scattering approximation. Namely, we need to include only contributions from processes 
shown in Fig^t], i.e. ladder diagrams, which amount to a self energy insertion [2| to the 
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photon propagator (JT5j) in the dominator as the following 

j,m,k' 

1 1, q; rrij] k, k') (1, q; rrij; k, k'| G^ - 1 (E) |1, q; m^; k, k') (1, q; rrij] k, k'| 
+ |1, q; rrij] 0) (1, q; rrij] 0| G (0) (22) |1, q; m 3 -; 0) (1, q; rrij] 0| 

xy|l,q;^;k') (18) 

Using the prescription — > N f <2R and Y^y ~* J^f I ^ we can cnan g e sums to 
integrals. The integration over Rj gives 5 (k — k') as a result of the translational invariance 
of the medium. The integration over k' picks up contributions only at k and the self energy 
can be written as 



+ 



(19) 



E — tku' k — ftw E — hw r k — hw — 2fko k 

Spontaneous decays of the medium atoms are ignored because medium atoms are far detuned 
from the source atom. The normalized photon propagator, with the above self energy 
modification, is 

- W = - W (20) 

E + i0 — fku r k — hwk — J2 E - hjjj r k - [I - \Na (E)] rvjj k ' ( J 

where the polarizability is defined as 

•<*>~*S (21) 



and the detunning 



= 1 (22) 

A (E) E — fvjj r k — flu E — tiw k — hw — 2hu k v 1 



The self energy insert ion (fT9*|) brings a correction of order Not to the denominator of the 
photon propagator, which cannot be obtained from any finite order calculation. Substituting 
this normalized photon propagator (J2*U|) back to the first order formula ()13|) and using the fact 
that W 2 = W, one fines the relevent reolvant matrix element to be 

(1, -hk; g m ; k; 0| G (E + zO) |2, 0; g m ; 0) 

1 1 

= %fc (/^S • e k) ^ Z ■ : r — — 7- 1 r-r (23) 

E — nw + innq E — hijj r k — (1 — %Na) fku k 
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To find the transition amplitude in Equation©, we need to integrate over E according to 
equation (J1UJ1. The integration includes contributions from two poles, one at E — fuv k + 
(l — \Not) fajJk, and one at E = hujo — ihnj. In the limit ^> 1, only the first pole 
contributes, since the second one has a finite imaginary part and its contribution decays 
away as e _nrr . The amplitude of finding a photon with wave number k and source atom 
with momentum — Kk is then 

7 ,. */ * exp [-i (hw r k + (1 - \Nat) u k ) r] 

h = hm g k (// • e k ) 1 r — ^ — — J - (24) 

ou k + (1 — 2^ a ) - + iwy 

and the corresponding probability is jl^] 

\bk\ 2 = \gk\ 2 |/Uo • £fc| 2 2 (25) 

[uJk/n - (wo - ii 2 uJq)] +n 2 7 2 



with cJq = JtlJq/2Mc 2 . We have identified n = 1 + \Na in the limit of small Na ll|. The 
photon frequency centers around n (u — n 2 f\WQ) . The total probability of finding a photon 
in the field is just the sum over all the k's 

y M 2 = ^lj^- r dUk ^ — — = i (26) 

V 3 (2vr) 2 Tic 3 7 K/n-^-nH)] 2 + n 2 7 2 

Note that in working out the integration, we have used the Wigner-Weisskopf approximation. 
We extend the integration of cu k to start from minus infinity |l4|, and take into account only 
the contribution from the pole at n (uq — u 2 uJq + iwy) 16]. Here the decay rate 7 should be 
evaluated at the photon frequency and it is given by 7 = 2fi 2 ujl/3hc 3 . 

Other relevent amplitudes, indicated in Fig[2], can be calculated by the same technique 
and now consider the processes shown in Fig[2]. Making use of the normalized photon 
propagator (|2*nj) . we find the amplitude for a medium atom located at Rj, being excited to a 
sublevel m, without any photon present to be 

lim T | gfc ]VV k ^ eXpH ; ( ^ +( !~^ a)a;fc)rl (27) 
- - u k + (1 - 2^") ^k - ^0 + inq 

-W m , (28) 



r— »oo ' J 

k 



(1 - \Na) UJ k - UJ 

It proves more useful to work in momentum space, since one can take advantage of the trans- 
lational invariance of the problem. We label the corresponding momentum state amplitude 
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FIG. 2: Medium atoms are excited by emitting or absorbing normalized photons. Figure (a) 
corresponds to a rotating situation while figure (b) corresponds to an anti-rotating situation. 



as b m (k) 



K <*) = U m « W » «p[-.- (M+(l-^aK)r] 

col + 1 - iiVo; )u k — Un + iwy 



X 



(1 - iJVa) uo k 
The probability of finding medium atoms excited is 



,, r (l - \Na) u k -L) Q + iwy 



\9k\ 



J2K(k)\ 2 = NV^J2r-/ ( 2 rtf+ 2 2 
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(29) 



(30) 



1 - \Na) cu k 

The amplitude of the anti-rotating excitation, a medium atom excited with two photons 
present can be shown to be 

exp [-i (Huj r k + (1 - \Na) uj k ) t] 



Knkk (k) = Jim VNV fi 2 \g k \ 2 

1 



t— >oo 



X 



(l - \Na) uo k -uo- 2uj, 
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(31) 
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and the corresponding probability, summing over different magnetic sublevels, 

i i 4 

17 11 \|2 ,r,r 4 \9k\ 



^\b mkk (k)\ 2 = Nv^r-, — t 

2 



[uJk/n — {u Q -ti 2 loI)\ +n 2 Y 

W ,o (32) 



(l — |iVa) ujk — uj — 2uJk 

The probability £) m ,ft |&m (&)| 2 and ^ mfc |& m fcA: (&)| 2 are of order Nau) k /uj, which are negli- 
gible with uo k <C u; and iVa -C 1. By now we have shown all the nonvanishing amplitudes of 
the system at r — ■> oo. The total probability, sum of equations (J2HJ) (JBOJ) (JH2J) , is found to be 
one to the order of Net. 

The average recoil energy of source atom is calculated solely from the amplitude b k ()24j) 

as 

This result shows that the source atom recoils according to the canonical photon momentum 
nko- This modification of the source atom recoil is directly related to the fact that the photon 
in the medium centers at a frequency nu Q + O {TidI/Mc 2 ) instead of u + O (/z^q/Mc 2 ), as 
is shown in formula (|24j) . In the spontaneous decay of the source atom, momentum of the 
source atom plus the field is conserved. This modification of the center freqency of the 
photon therefore results in a modification of the source atom recoil. However, this shift 
in the photon frequency seems strange, because we expect the frequency to center around 
the atomic frequency Uq from a energy conservation consideration. In the following, we 
give a detailed analysis of the energy conservation to order of the shifted frequency, namely, 
nhw — hujQ = \NahjjQ. 

Before we proceed to show energy conservation, it is necessary to note that the source 
atom excited at t — is not an eigenstate of the system. However, a discussion of energy 
conservation is meaningful in the context of average energy being conserved. 

In order to find all the other forms of energy besides the photonic excitation, we should 
note that as the photon moves along in the medium, medium atoms inside the sphere of 
R = ct are excited, according to e~ n/ ^ ct ~ R \ Even the excitation probability is of order of 
magnitude Naujk/oj, as far as the energy is concerned, this produces a correction of order 
NauJu- Including this contribution to the energy, and the interaction energy of the medium 
atoms with the field, we should be able to recover the conservation of energy. In the following 



11 



we work out explicitly energies associated with different excitations in the system and show 
that the total average energy is conserved. 

The first part of the energy corresponding to the field excitation is 

1 



y] \h\ 2 f^k 



1 + -Na ) hw 



(34) 



The one associated with the medium atoms excitations is 



^ \b m k\ 2 + \bmkk\ 2 (hio + 2hu> k 
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— — = — hui Ti 
A' 



1 



-1 2 



+ 



1 



(l + iJVa) hw k + huj 



(l — \Na) hu k — hui 

Here the difference in the energies of state b mk and b m kk can be ignored, because it amounts 
to a correction in the order Nau k /uj. To leading order of uik/uim, we find that 

1 



1 
A 7 



A(E) 



—^—Na— 

hw uj 



(36) 



which enables us to neglect the difference between and ^p^y- Here the detuning A (E) 
is evaluated at the pole of the photon propagator E — (l — \Na) hio k . Making use of the 
definition (|21jl . one finds that the energy associated with the atomic excitations is ^NahujQ. 
The third part of the energy, the interaction energy between excited medium atoms and the 
field can be calculated as, 



(V) = N V h 9k (Kbm + Kbmkk) + h.c. 



m,k 



} J -Na (ui k ) hui k - 



\hg k 



(hwk/n — hu>o + n 2 hu> r k ) + n 2 h 2 ^ 2 
-Na (u> k ) hj k 



(37) 



adding up all of these contributions, we have 

^2 \bk\ 2 hui k + ^2 (\ b mk\ 2 + \b m kk\ 2 ) hu+ (V) = huj (38) 

k m,k 

showing that the frequency shift of the photon in the medium is compensated by the exci- 
tation energy of medium atoms and the interaction energy of excited atoms with the field. 
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III. DISCUSSION 



In this paper we have shown that the source atom recoils according to nhko, which agrees 
with the previous theoretical and experimental results This modification of the 

photon recoil arises in our calculation as a result of the interference of the different scattering 
amplitudes of the photon. As has been shown in the calculation, while the center frequency of 
the soure atom is at u> , only frequencies centered at nu experience contructive interference. 
This is very much similar to the case of a source atom radiating in a cavity with the cavity 
frequency instead of the atomic frequency. From a quantum point of view, the source 
atom decays because it radiates and reabsorbs virtual photons. Such a process introduces 
a finite self energy whose real part give s the level shift and whose imaginary part gives the 
decay of the atomic excitation |?| This process, though not included explicitly in our 

calculation, is the only way that the source atom "knows about" the enviroment (vacuum, 
cavity or dielectric medium) in which it locates. In the cavity, the virtual photon can 
interfere with those reflected from the cavity walls and constructive interference occurs only 
at the cavity freqency. In the dielectric medium the virtual photons radiated by the source 
atom can be scattered by medium atoms and reabsorbed by the source atom during the 
time t z< I/7. Different scattering amplitudes interfere to shift the real radiating frequency 
to nu . This is a different effect from the level shift, the real part of the source atom self 
energy, due to the interaction of the source atom with the enviroment. In our particular 
example of the source atom located in a dilute medium, the shift due to interacting is of the 
order Nof) while the shift due to interference is of the order NaujQ. 

An alternative explaination can be put forward in terms of eigen excitations of the 
system. Actually, if we consider the interaction of the field with the medium atom, the 
eigenmode of the system is neither medium atoms being excited nor a photon present, but 
a superposition of these two type of excitations, or a polariton Q]. The source atom decays 
by radiating polaritons instead of photons. For the energy to be conserved, the polariton 
energy plus the recoil energy should equal to the initial average energy hujQ. As we have 
shown, the photon carries only part of the energy of the eigen excitation. On the other 
hand, the momentum carried by the medium atoms is negligible, the polariton momentum 
is just the photon momentum. When we require energy and momentum conservation for 
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the radiating process, the only possibility is the source atom to recoil according to nhk . 
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